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Abstract: In the 3- flavour framework we derive a simple approximate analytic expression 
for the day-night difference of the flux of solar z/g at terrestrial detectors which is valid for 
an arbitrary Earth density profile. Our formula has the accuracy of a few per cent and 
reproduces all the known analytic expressions for the Earth matter effects on the solar neu- 
trino oscillations obtained under simplifying assumptions about the Earth's density profile 
(matter of constant density, 3 layers of constant densities, and adiabatic approximation). 
It can also be used for studying the Earth matter effects on the oscillations of supernova 
neutrinos. We also discuss the possibility of probing the leptonic mixing angle ^13 through 
day~night asymmetry measurements at future water Cherenkov solar neutrino detectors. 
We show that, depending on the measured value of the asymmetry, the current upper 
bound on ^13 may be improved, or even a lower bound on this mixing parameter may be 
obtained. 
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1. Introduction 

Solar neutrinos coming to terrestrial detectors at night travel some distances inside the 
Earth and so their oscillations are affected by the Earth's matter. This leads to a difference 
between the day-time and night-time solar neutrino signals - the day- night, or "regenera- 
tion" effect [1-7]. Solar neutrino oscillations ^ depend mainly on two parameters, the mixing 
angle 612 and the mass squared difference Amg^. For example, the recent three-neutrino 
global fit of solar, atmospheric, reactor and accelerator neutrino data of Ref. [9] gives for 
the solar neutrino oscillation parameters the 3a allowed ranges 9i2 = (28.6 -7- 38.6)° and 
Am|i = (5.4^9.5) x 10"^ eV^, with the best-fit values Ou = 33.2° and Aml^ = 6.9 x IQ-^ 
eV^ . In the three-fiavour framework there is an additional dependence on the mixing angle 
^13 which determines the component of the electron neutrino in the third mass eigenstate 

^The presence of non-oscillation phenomena is now strongly constrained by a combination of KamLAND 
and solar neutrino data, see, e. g. Ref. [8] and references therein. 
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ly^ separated by a large mass gap from the other two mass eigenstates. This mixing pa- 
rameter is mainly constrained by the data of the CHOOZ reactor experiment [10] and is 
known to be small, 

^13 < 9.8° (13.4)° , or sin ^13 < 0.17 (0.23) (1.1) 

at 90% C.L. (3cj) [9]. The smallness of ^13 implies that its effect on the solar neutrino 
oscillations is rather mild; however, with increasing accuracy of the data even relatively 
small effects have to be taken into account. 

Three-flavour analyses of the solar neutrino data, including the day- night effect, have 
been performed in a number of recent studies [9, 11-15]. In particular, in Ref. [15] a 
simplified analytic treatment of the day-night effect in the three-flavour framework was 
given, and it was shown that the day-night z/g flux difference scales essentially as cos® ^13 
rather than cos '^^13, which is the expected scaling law for the day-time flux. The authors 
of [15] considered the Earth matter effect on the solar neutrino oscillations under the 
assumption of constant matter density of the Earth. The purpose of the present paper is 
to extend their analysis beyond the constant density approximation, and to scrutinize to 
what extent day-night asymmetry measurements may become a significant way to probe 
for the mixing angle ^13. Using the relative smallness of the matter-induced potential of 
neutrinos inside the Earth, we derive a simple and accurate analytic formula for the day- 
night flux difference in the case of an arbitrary Earth's density profile. We compare our 
result with the known analytic formulas obtained under simplifying assumptions about the 
density profile of the Earth as well as with the results of the exact numerical calculations. 

The paper is organized as follows. In sec. ^ a general consideration of the day- 
night effect in the 3-flavour framework is given. In sec. ^ we present our main result 
- the derivation of an approximate 3-flavour analytic formula for the day-night i^e flux 
difference in the case of an arbitrary Earth density proflle. We also discuss a number 
of special cases for which the day-night flux difference can be found analytically. These 
include approximating the matter density profile of the Earth by one layer or three layers 
of constant densities, and also the adiabatic approximation to neutrino evolution inside the 
Earth. In sec. |^ the comparison of our formulas with the results of the exact numerical 
calculations is performed. In sec. ^ we make a quantitative study of the simple correlation 
between the day-night asymmetry in the solar neutrino flux and the magnitude of the 
mixing angle ^13 characterizing leptonic CP violation in neutrino oscillations ^. We use 

^ There are additional CP violating phases but they do not afltect neutrino oscillations, only lepton- 
number changing processes, like neutrino-less double beta decay [16-19]. 
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this correlation in order to forecast the degree to which the mixing angle ^13 can be probed 
at future solar neutrino measurements of the day-night asymmetry. We discuss our results 
in sec. |6|. In the Appendix, we present the derivation of an improved perturbation-theoretic 
formula for the Earth regeneration factor. 

2. Day-night effect with three flavours: Generahties 

The 3-flavour day-time probability of finding an electron neutrino produced in the Sun in 
the electron neutrino state at the detector is to a very good accuracy given by [20] 

Pd = 4^^2/(^12, Amii, 4y) + 4 . (2.1) 

Here P2f{0i2, Am^i, cf^V) is the 2-flavour survival probability in the Sun, calculated 
with the neutrino potential in the Sun V substituted by the effective potential cf^V, and 
we use the standard notation Sij = sin 6ij , Cij = cos 9ij . Upon averaging over fast neutrino 
oscillations inside the Sun, this probability can be written as 



1 



P2f{ei2,Aml^,cl^V) = - [1 + cos 2012 cos 2012 (1-2P') 



(2.2) 



where P' is the probability of transition between the first and second matter eigenstates 
in the course of the neutrino propagation in the Sun ("hopping" probability), Ouir) is the 
effective mixing angle in matter at the point r, and cos 29i2 is the value of cos 20i2 averaged 
over the neutrino production point [15]: 



Q 



cos 26*12 =/ cir/(r) cos 2012 (r) . (2.3) 







Here /(r) is the normalized spatial distribution function of the neutrino source in the 
Sun [21], and the explicit formula for cos20i2(r) will be given below [Eqs. ( p.l6| ) and 



The solar neutrino fiux arriving at the Earth is an incoherent sum of the fiuxes of 
mass-eigenstate neutrinos vi, V2 and 1^3 (see, e.g., [22]); therefore the day-time V(, survival 
probability can in general be written as 

Pd = PH PfJ + Pf2 pS^ + P^s P£^ , (2.4) 

where P^ {i = 1,2,3) is the probability of u^. — > i/j conversion in the Sun, and pj^^^ is the 
projection of the ith mass eigenstate onto z/g: P^^^ = {Ueil"^, U being the leptonic mixing 
matrix in vacuum. For this matrix we use the standard parameterization [23] 

u = 023rsOi3rlOu 
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-S12C23 



C12C13 

- Cl2Sl3S23e"^<^P C12C23 



S12C13 



5136 



-i<5f 



4<5, 



V S12S23 - Cl2Sl3C23e 



iSr 



-C12S23 - Sl2Sl3C23e 



(2.5) 



C13C23 / 



where Oij is the orthogonal rotation matrix in the ij-plane which depends on the mixing 
angle 9ij, and Ts = diag(l, l,e^'^cp), 6cp being the Dirac-type CP-violating phase. 

As we shall see shortly (see also the discussion in Ref. [15]), the third matter eigenstate 
u^rn inside both the Sun and the Earth essentially decouples from and f2m- In addition, 
it is practically not affected by the solar or Earth matter, i.e. — 1^3- Therefore one 
5^3. It is also easy to show that upon averaging over the fast neutrino 



has P^^ ~ R 



,{0) 



3e 



oscillations in the Sun one obtains [15] 







'13 



l + cos26'i2(l-2P') 



p0 



rl3 
2 L 



l-cos 26*12 (1-2P') 



(2.6) 



Substituting ( |2.6D into (^^ , one recovers Eq. ( ^ ) with P2 / (6*12 , Am2i , cf 3 V) given by 
Eq. ( |2.2D . Note that for the LMA MSW solution of the solar neutrino problem P' is 
extremely small; therefore we shall neglect it from now on. 

Since the flux of solar neutrinos coming to the Earth is an incoherent sum of the fluxes 
of mass-eigenstate neutrinos, the night-time survival probability at the detector can be 
written as 

= Pi Pt + Pi P^e + P% Pg , (2.7) 

where P-® is the probability that a neutrino arriving at the Earth as a mass eigenstate 
is found at the detector in the V(, state after having traveled a distance L inside the Earth. 
The night-day probability difference is therefore 



P, 



N 



Pd = P^i {P?e - P?^) + P% {Pi - P^h + P% {Pie - Plt^' 



The probabilities P^, P-^ and P^^^ satisfy the conditions 



i=l 



i=l 



i=l 



which follow from the unitarity of neutrino evolution. Since P 



3e 



4^ Eq. 



pfcB pi. 



(0) 



= -(P2^ - P^^^). Substituting this into Eq. (|2j) yields 

Pn-Pd = {Pi - P!i){Ple - P^'J) = -C?3COs2^i2(P2®e " P 



(Oh 



2e J ' 



(2.8) 

(2.9) 
I) gives 

(2.10) 



where in the last equality Eq. (2.6) with P' = has been used. The Earth matter effects 

2e 



on the solar neutrino oscillations are encoded in the factor P® — P2g^ ; we now turn to the 



calculation of this quantity. 
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Neutrino evolution equation in the flavour basis can be written as 

f/diag(0, 2(5, 2A) U"^ + diag{V{t), 0, 0)1 i' , 



. d 
i—u 
dt 



(2.11) 



where v = {ve, f^, Vt)'^ is the wave function of the neutrino system, t is the coordinate 



along the neutrino trajectory, and the parameters 

6 = ^, A 



(2.12) 



AE ' AE 
involve the mass-squared differences responsible for the solar and atmospheric neutrino 
oscillations, with Amgj^ large compared to Arn^i [9]: Amii ~ 2 x lO'^ eV^ 
The matter-induced potential of neutrinos in Eq. ( |2.11D is 



V{t) = V2GFNe{t), (2.13) 

where Gp is the Fermi constant and N^^t) is the electron number density in matter. 
For our purposes it proves convenient to go to the new basis defined through 



The evolution equation for the rotated P state is 

(isj^S + cj^Vit) 2S12C126 Si3Cl3V{t) \ 



(2.14) 



. d 

I— 
dt 



(2.15) 



2si2Ci2 5 

\ Sl3Cl3V{t) 2/^ + V{t)sl^J 

For neutrino evolution in the Sun and in the Earth one has V ^25 2A (actually, inside 
the Earth V <C 25). Since, in addition, S13 <C 1, one can to a very good accuracy neglect 
the (1-3) and (3-1) elements of the effective Hamiltonian in Eq. ( |2.15| ) compared to the 
(3-3) element. This means that the evolution of the third matter eigenstate essentially 
decouples from that of the first two eigenstates and, in addition, that matter effects on the 
third eigenstate are negligible. This approximation is especially good for the survival 
probability since, as was shown in [24, 25], terms of first order in S13 cancel in this quantity, 
and the correction only appears at order Si^V"^ / A"^. 

Diagonalization of the (Pi, D2) - subsector of the effective Hamiltonian in Eq. ( 2.15| ) 
yields the instantaneous effective mixing angle in matter ^12 (i), which can be determined 
through 



cos 26112 (t) 



cos 2012 ^ 



-13 



V{t)/2 



uj{t) 



where 



co{t) = J [cos 2^12 6 - cf3y(t)/2]2 + 52 gin2 2612 . 



(2.16) 
(2.17) 
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Let us introduce the neutrino evolution matrix in the rotated basis according to 



i>(t) = S{t, to)i){to) , 



Sito, to) = l. 



(2.18) 



The matrix S{t, to) satisfies the same evolution equation, Eq. (|2.15| ), as u. The decoupling 
of the third eigenstate implies that it can be written as 

/ a{t,to) Pit, to) \ 



Sit, to) 



\ 



(2.19) 



-P* it, to) a* it, to) 
fit,to)J 

where 

/(t,to) =exp[-2iA(t-to)] , (2.20) 

and the parameters a(t, to) and /3(t, to) (satisfying |ap + = 1) are to be found from the 
solution to the two- flavor problem governed by the mass squared difference Am^^, mixing 
angle 9i2, and matter potential cl^Vit). 

It is now easy to find the probabilities and in terms of a and (3. Direct 
calculation yields P^® = c1^\ci2C( — Si2/?P, P^ = cf3|si2a + ci2/3p, and 



K - He = c?3 [cos 2^12 \(3\'' + sin 2^12 Re(d*/?)] . 



(2.21) 



3. An analytic expression for P^r — Pd 



We shall now derive an approximate analytic expression for P^ — Pg^ ^ valid for an arbitrary 
matter density profile of the Earth. Our basic point is that the neutrino potential in the 
Earth V is small iV/26 < 0.05), and so can be considered as a perturbation. 

We shall use the standard perturbation theory for the evolution matrix (see, e.g. Ap- 
pendix B of Ref. [25]). The effective Hamiltonian in the rotated basis in Eq. ( ^.15 ) can be 
decomposed as Hit) = Hq + Hiit), where Ho and ^i(t) are of zeroth and first order in 
Vit), respectively. Then to first order in Vit), the evolution matrix Sit, to) can be written 
as ^ 

Sit, to) Soit, to)-iSoit, to) [ [Soit' ,to)-^Hiit')Soit' ,to)]dt' . (3.1) 

J to 

Both Soit, to) and 5(t, to) have the form ( 2.ig| ). For the zeroth-order matrix S'o(t, to) we 
find 

ao = cos 5L + i cos 26*12 sin 5L , Po = —i sin 26i2 sin 5L . (3.2) 
Eq. (|3.1| ) then gives 



Q = ao(l — iA) — iPoB* 



(3 = j3oil + iA) -iaoB, 



(3.3) 
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where 

^ = ^ / (l«oP - |/3oP) Vdt' , B = clJ' d*Jo Vdt' . (3.4) 

Substituting Eqs. (U) and (U) into Eq. (|2^ ), one arrives, after a little algebra, at a 
remarkably simple result: 

1 

P®e - P2e = Ci3 sin^ 2^12 - dx V{x) s\n[25 ■ (L - x)] . (3.5) 

2 Jo 

Here x is the coordinate along the neutrino path inside the Earth and L = 2Rqcos6z, 
where is the Earth's radius and 9z is the zenith angle of neutrino trajectory. The 
night-day flux difference Pn-d is then given by (|2l0| ). Equations (|2lo|) and (U) are our 
main results. 

Strictly speaking, for perturbation theory in V to be valid, two dimensionless param- 
eters have to be small, namely V/25 and VL (or, more precisely, Vdx). While the first 
parameter is indeed very small in our case, for long enough distances travelled by neutrinos 
in the Earth the second parameter may be rather large. Nevertheless, as we shall see in 
sec. Eq. ( |3.5D yields very accurate results even for large VL provided that integration 
over neutrino energies or zenith angles is involved. A more accurate formula can be ob- 
tained by replacing in the integrand of Eq. ( |3.5| ) the in-vacuum oscillation phase by the 
corresponding adiabatic phase: 

P® - P!f^ = sin^ 2^12^ ^ dx V{x) sin 2 j w(x') dx' 



(3.6) 



where uj{x) is given in Eq. ( pTl) . We shall show in sec. H that this considerably improves 
the agreement with the exact results in the cases when the oscillation phase is large but 
no averaging over zenith angles or neutrino energies is performed. Eq. ( |3.6D is based on 
perturbation theory in V/25 rather than in V and thus is valid for arbitrary values of VL. 
We give its derivation in the Appendix. 

It is interesting to note that Eq. (|3.6| ) depends on the adiabatic phase even though no 
explicit assumptions about adiabaticity is made in its derivation. This is related to the fact 
that in the limit V -^25 the adiabaticity condition is satisfied automatically, irrespective 
of the spatial behaviour of the potential V{x). Indeed, the adiabaticity condition requires 
that the mixing angle in matter 612 change little on the length scales of the order of the 
neutrino oscillation length in matter. In the case V <^25 this mixing angle is always close 
to the vacuum mixing angle 612 and so changes little over the whole scale L which contains 
many oscillation lengths, i.e. the neutrino evolution is automatically adiabatic. 
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We will consider now several special cases for which exact analytic expressions can be 
found for a and /3 (and so also for P®^ — P^)- 

The matter density of the Earth changes relatively slowly both within its mantle and 
core, while at the mantle-core border it jumps by about a factor of two. Therefore as a 
first approximation one can consider the matter densities in the mantle and core as being 
constant. In other words, for neutrinos traversing only the mantle of the Earth one can use 
the constant density approximation, whereas for core-crossing neutrinos the Earth density 
can be approximated by three layers of constant densities (mantle-core-mantle). More 
accurate treatment along these lines would include a larger number of layers, allowing for 
several layers within the mantle and the core [26]. 

For existing solar neutrino detectors, most of the neutrinos detected during the night 
time cross only the mantle of the Earth but not its core (the Super-Kamiokande (Super-K) 
detector has the largest fractional core coverage time equal to 7%). We therefore start with 
the case of one layer of constant density. 

3.1 Matter of constant density 

In this case one readily finds 

a = cosojL + icos29i2 smojL , /? = —i sin 2^12 sin , (3.7) 

where the mixing angle in matter 9i2 and lo are given by Eqs. (|2l|) and (gl^) with 
y=const. Substituting (3/7) into Eq. (|2.2l| ) and the result into Eq. (^.10|) , one finds 



P^e - P2e = 4 sm' 2012 ^ siu^ , (3.8) 

V5 



Pj^- P^ = -c1. COS 2012 sin^ 26*12 — k sin^ ujL . (3.9) 

This is the result obtained in Ref. [15]. In the limit S13 one recovers the 2-flavour day- 
night probability difference in the constant-density approximation to the Earth's density 
profile given in [27]. To leading order in the potential V Eq. ( |3.9D simplifies to 

V 



Pj^- P^ = -c^g cos 2^12 sin^ 26*12 77? sin^ 5L . (3.10) 
3.2 Three layers of constant densities 

For neutrinos crossing the core of the Earth, the Earth's density profile can be approximated 
by three layers of constant densities pi, p2 and pi and widths Li, L2 and Li (mantle-core- 
mantle). In this case one can write a = Z — iW^, (3 = —iWi with real Z and VFi,3, which 
yields P^^ - P^°^ = cj^Wiicos 2612 Wi + sin 2^12 W^). Here [28] 

Wi = 2 sin 2^1 Y sin coiLi + sin 2^2 sin LO2L2 , 
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Ws = —(2 sin 2^1 Y cosuJiLi + sin 2^2 cosw2-^^2 ) , 



(3.11) 



with 

Y = coswiLi cosa;2-Z>2 — cos(20i — 2^2) sinwiLi sinci;2i2 , (3-12) 

and 61 and 62 are the values of the effective mixing angle 612 in matter of densities pi and 
P2, respectively. After a simplification one then obtains 

^2® - ^2?^ = 4 sin2 2^12 — (2y sinc^iLi + — smuj2L2) 

ZlJ\ U!2 

X ( — 2y sinwiLi + — sina;2i^2 ) • (3.13) 
V^i L02 ) 

In the 2-flavour limit ^13 — > the expression for P®^ — for the 3-layer model of the 
density profile of the Earth was first obtained in [29] . In the limiting cases V2 = Vi, L2 = 



or Li = Eq. ( 3.13 ) goes into the corresponding expressions for one layer of constant 



density, Eq. (3.8). 



To first order in Vi^2 Eq. ( [3.13 ) simplifies to 



P2e-P2e = 4 sin^ 26^2 ^ [Vi sin5(2Li + L2) + {V2 - Vi) sin5L2] sin5(2Li+L2) . (3.14) 
3.3 Adiabatic approximation 

If the matter density changes slowly enough along the neutrino path, the evolution of the 
neutrino system can be studied in the adiabatic approximation. This approximation should 
be especially good for neutrinos crossing the Earth's mantle only, and can be considered 
as an improvement of the constant-density approximation. For core-crossing neutrinos one 
can consider the evolution inside the core and inside the mantle as adiabatic and match 
the wave function of the neutrino state at the mantle-core borders; here we discuss the 
adiabatic evolution for mantle-crossing neutrinos. 
In the adiabatic approximation one has [24] 

a = cos{6f — 6i) coscf) + icos{6f + 6i) smcp , 

P = sm{6f — 9i) cos</> — isin(6'/ + 9i) smcf) , (3.15) 



where 



uj{x)dx , (3.16) 







and 6i and Of are the values of the mixing angle in matter 9i2 at the initial and final 
points of neutrino evolution. To a good accuracy, the matter density profile of the Earth 
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is spherically symmetric, so that one can set Of = Of, however, for generality we shall allow 
for different values of 9i and Of. 



Substituting ( 3.15 ) into Eq. ( p.21| ), one finds 
P2e-P2e = 4 [^H^ f " ^i) f - Oi + 2^12) COS^ + sin(^j + Oi) sin(^/ + 0, 



To leading order in Vij this simplifies to 



2e 



P2^^=cf, sin2 20 



13 



yi2 



cos oL H sm oL 



25 



25 



20i2)sin2(/)] 
(3.17) 

(3.18) 



It is easy to make sure that for all the special cases discussed above our formula ( |3.5| ) 
exactly reproduces the corresponding leading-order in V expressions. The cases of one 
layer of constant density and three layers of constant densities [Eqs. ( |3.10| ) and ( 3.1^ )] 
are immediately obtained from Eq. ( p.sp by a straightforward calculation of the integral 
with the corresponding matter density profiles. The adiabatic case is recovered when one 
integrates in Eq. (^^) by parts and neglects the term with the derivative of the potential 
(i.e. keeps only the off-integral term). This gives 



2e 



^2e — 



1 



4, Sin' 20u ^ {Vir) cos[2J • (L - r)]} 
= ^3 sin^ 26*12 {^f-^i cos 25 L} , 



(3.19) 



which coincides with Eq. ( p. 18 ). 

Quite analogously, in all these special cases Eq. ( |3.6|) exactly reproduces the corre- 
sponding leading-order in V/25 expressions for — ■ 



4. Comparison with numerical results 

We now proceed to test the accuracy of our approximations by comparing the results 
obtained using our main formulas in Eqs. ( 2.1[1| ) and ( |3.5| ), as well as Eq. (p.6|), with those 
following from the exact numerical integration of the three-neutrino evolution equation 
employed in Ref. [9]. 

As has already been pointed out, for the validity of perturbation theory in V the con- 
dition V -^25 \s not sufficient: one must also require that the correction to the oscillation 
phase, which is ~ VL, is small compared to unity, and not only compared to the main 
contribution to the phase, 5L. This is illustrated in Fig ||, which gives the comparison of 
the predicted probability differences due to the regeneration effect at the Earth, P®^ ~ P2e ^ 
versus zenith angle, corresponding to different solar neutrino arrival directions. The solid 
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0.08 




Figure 1: Two-flavour Earth regeneration factor versus zenith angle for a 10 MeV solar neutrino 



from the exact numerical calculation and the analytic form given in Eq. (3.5) 



(red) curve follows from the exact numerical calculation, while the dashed (black) curve 
corresponds to the analytic formula in Eq. ( |3.5| ). One sees clearly that the difference be- 
tween the night-day probability differences found in our analytic approximation and those 
that follow from the exact numerical integration of the three-neutrino evolution equation 
is maximal for neutrinos with small zenith angle. For such neutrinos with a longer path 
inside the Earth the difference between the exact and approximate oscillation phases accu- 
mulates and becomes relatively large, leading to significant differences between the exact 
and analytic results, as can be seen from the figure. 



0.04 



0.03 



-0.01 



numerical 
analytical 




0.2 



0.4 0.6 
COS 6 



0.8 



Figure 2: Two-flavour Earth regeneration factor versus zenith angle averaged for the Super-K 
detector, from the exact numerical calculation and the analytic formula in Eq. (|3.5|). 
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On the other hand, for large enough baselines, when VL becomes of order unity, also 
6L ^ 1, so that the averaging regime sets in, either because of the integration over the 
baselines or, for a fixed zenith angle, because of the integration over energies. Note that 
averaging over energies is always present in the experimental data due to the fact that the 
neutrino energy is not directly measured and also due to the finite energy resolution of the 
detectors. In addition to the energy resolution, for comparison with data one has also to 
fold in the neutrino spectrum and relevant neutrino-detection cross section. All this leads 
to the smoothening of the above difference. The net effect is that even order 1 shifts of the 
phase become unimportant in practice, and neglecting the matter-induced correction to the 
phase is justified. This is illustrated in Fig. ^, which compares the probability differences 

— plf^ properly averaged over the full energy range relevant for the Super-K detector. 
It can be seen from this figure that the averaging gives a much better agreement between 
the analytic and exact results. 

0.08 



0.06 



0.04 - 



Oh 0.02 



-0.02 



numerical calculation 
adiabatic improvement 




0.4 0.6 

COS 9 



Figure 3: Two-flavour Earth regeneration factor versus zenith angle for a 10 MeV solar neutrino 



from the exact numerical calculation and the improved analytic formula given in Eq. (3.6) 



Finally, note that due to relatively low statistics in the solar neutrino experiments, 
integration over zenith angles is also usually required. This also leads to an averaging 
over the oscillation phase, giving a very good agreement between the analytic and exact 
numerical results (see Fig. ^ below). 

In the situations when no (or little) averaging over neutrino energies or integration 
over zenith angles is involved, the improved perturbation theory formula given in Eq. ( |3.6| ) 
yields much better results. This is illustrated by Fig. |^ which shows that in this case there 
is perfect agreement between the analytic and numerical results. 
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Figure 4: Comparison of predicted Earth regeneration factors and day-night survival probabihty 
differences: sohd lines correspond to exact numerical calculation, dashed ones are the results of the 



analytic formulas in Eqs. (2.1C) and (3.5). See text for details. 



Finally, we turn to the case of interest of three-flavour day-night solar flux differences. 
Fig. ^ compares, for the site of the Super-K detector, the analytic and numerical predictions 
for the probability differences 

P2® - -Pi? and Pn-Pd, integrated over the whole night and 
day ranges of the zenith angles and averaged over one year. The results are given for the 
neutrino energy E = MeV. 

The left panels give the dependence of these probability differences on the angle ^13. 
The middle and right panels display the dependence on the "solar" mixing angle 612 and 
solar mass-squared difference Am^]^, respectively. The vertical bands in the left and middle 
panels indicate the allowed 3ct regions for sin^ ^13 and sin^ 612 respectively. The central 
curves in all panels are obtained for the best-fit values of the undisplayed neutrino os- 
cillation parameters, as obtained in Ref. [9]. On the other hand the outer curves in all 
panels correspond to the currently allowed Scr range of the corresponding undisplayed pa- 
rameter. For example in the left panels the upper curves correspond to sin^ 612 = 0.39 
and Am^ = 5.4 x 10~^ eV^, while the lower ones correspond to sin^ ^12 = 0.23 and 
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Am^ = 9.5 X 10"^ eV^ for energy E = 10 MeV. In other words, the separation band 
between the upper and lowermost curves in each panel gives a measure of the uncertainty 
in the theoretical prediction allowed by current data. 

Clearly, in all cases our approximate results reproduce very well the results found from 
the exact numerical integration of the three-neutrino evolution equation, with a precision 
much better than current experimental sensitivities. 



5. Future experiments and sensitivity to 9 



13 



The three-flavour day-night asymmetry in the solar neutrino flux correlates with the mag- 
nitude of the leptonic mixing angle ^13 [15]. Adopting the standard definition 

N -D 



N + D 



(5.1) 



and using ( |2.10| ) and (^) one can see that the main dependence on cos 6*13 is quadratic. 
Here we give a quantitative study of this interesting correlation, as displayed in Figs. ^ and 
|6[ The left panel in Fig. ^ gives the magnitude of the day-night asymmetry measured at the 
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Figure 5: Day-night asymmetries in three- flavour solar neutrino oscillations versus the mixing 
angle ^13. The left panel shows the present status, while the right one gives the projected sensitivities 
in future water Cherenkov experiments assuming the same central value of And. 

Super-K experiment 2.1% it 2.0%(stat) it 1.3%(syst) [30] compared with the range for this 
quantity which is theoretically predicted on the basis of current analysis of solar neutrino 
data [9]. The horizontal band corresponds to the current la uncertainty in the measured 
day-night asymmetry, while the hatched region shows the 3a uncertainty implied by the 
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current global determination of oscillation parameters. Clearly, the large errors with which 
the small day-night asymmetries are determined by current experiments are such that the 
asymmetry measurement does not play a significant role in constraining the angle ^13, the 
upper bound on which is dominated by other data in the global fit, mainly by the reactor 
data. However, the simple and direct correlation between the predicted asymmetries and 
^13 that can be appreciated in Figs. |5| and ^ can be used to provide an estimate of the degree 
to which the mixing angle ^13 can be probed at a future generation of water Cherenkov 
solar neutrino experiments. The significance of the future data for the determination of ^13 
will depend substantially on both the central value and the errors of the measured A^r^). 
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Figure 6: Same as the right panel of Fig. but for different central values of expected measured 
And- 2.8% (left panel) and 1.4% (right panel). 



The horizontal bands in the right panel of Fig. ^ and in both panels of Fig. ^ pro- 
vide an estimate of the expected improved day-night asymmetry measurements at a larger 
Super-K-like water Cherenkov detector, with la errors at 25% and 10% of the current 
Super-K error. Such errors might be realistic at future detectors. For example, the Un- 
derground Nucleon decay and Neutrino Observatory (UNO) experiment is expected to 
improve the sensitivity on the day-night asymmetry down to cr = 0.25 asK, while for the 
Hyper-Kamiokande experiment the errors can be much smaller [31]. 

On the other hand, the hatched regions in the right panel of Fig. ^ and in Fig. ^ cor- 
respond to the 3a range for solar neutrino oscillation parameters expected after improved 
KamLAND results. Although further improvements on the determination of solar neutrino 
oscillation parameters are likely to be available by then, the figure shows how the uncer- 
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tainty in the day-night asymmetry measurement dominates the projected accuracy on the 
mixing angle ^13. 

Let us first assume that the central value A]\fo measured in future experiments 
coincides with the current value found by Super-K (right panel of Fig. ^) . As one can see 
from the figure, no improvement on the current upper bound on ^13 can be obtained in 
that case. However, it is quite possible that the future central value of the measured night- 
day asymmetry will differ from the currently measured one, at least within the current 
la experimental error ^. This situation is illustrated in Fig. |6[ As one can see from the 
left panel, if the future central value of measured ^jv/) is higher than the present Super-K 
one, the current upper limit on ^13 can actually be improved. If, on the contrary, a lower 
value of And is measured (right panel of Fig. ^), the derived upper limit on ^13 will be 
substantially weaker than the current one and thus irrelevant. However, as can be seen 
from the figure, in that case a lower bound on ^13 may appear; together with the current 
upper bound it may actually lead to a rather precise determination of ^13. 

Thus, an improved measurement of the day-night asymmetry can be immediately 
converted into more precise information on ^13. 

It should be noted that, even though the day-time solar neutrino signal depends on 
613 more sensitively than the day- night asymmetry (essentially as cf^), this dependence is 
always multiplied by the overall normalization of the solar neutrino flux. For neutrinos, 
the uncertainty in the flux normalization factor fs is of the same size as the possible effect 
of non-zero ^13 . This makes it difficult to disentangle the two effects and hinders the precise 
determination of both fs and ^13 from the charged-current data. In contrast to this, the 
day-night asymmetry is independent of the overall flux normalization and so may provide 
an unambiguous information on ^13. 

6. Discussion and outlook 

We have derived a simple and accurate analytic expression for the day-night difference of 
the flux of solar coming to terrestrial detectors, with 3-flavour effects fully taken into 
account. Our approach was based on a simple perturbation theory in the matter- induced 
potential of neutrinos inside the Earth, without any assumptions regarding the Earth's 
density profile. Our results are therefore valid for an arbitrary density profile. 

^Future best-fit values of tfie parameters 612 and Ami2, which determine the position of the hatched 
regions m Figs. I and §, can also differ from the current ones. However, the effect of this deviation cannot 
be as large as that for And since their present-day errors are significantly smaller. 
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We have checked our analytic formula by comparing it with the results of the exact 
numerical calculations with the 3- flavour evolution equation, and found that the accuracy of 
the analytic approach is typically about a few per cent when the neutrino path length inside 
the Earth is small, or when the integration over the zenith angles or averaging over neutrino 
energies are performed. On the other hand, when both the neutrino energy and zenith angle 
are fixed, the formula in Eq. ( |3.5D does not provide a sufficient accuracy. For this case we 
have derived an improved perturbation-theoretic formula [Eq. (|3.6D]. Compared to the 



expression in Eq. (3^) it contains the adiabatic oscillation phase in the integrand instead 
of the in-vacuum one. The improved formula gives a perfect agreement with the exact 
numerical results even when no averaging over the zenith angles or neutrino energy (i.e., 
over the oscillation phase) is performed. 

We have studied the dependence of our results on the leptonic mixing angle ^13 which 
is of great interest for many reasons, most importantly because it governs CP violation in 
neutrino oscillations. We have found that for an arbitrary Earth density profile the day- 
night difference of the solar neutrino flux at the detector scales mainly as cf^ , as previously 
found for the constant density profile in [15]. The remaining (milder) dependence stems 
primarily from the 0i3-dependence of the neutrino mixing angle in matter at the production 
point inside the Sun. 

Although the smallness of 613 implies that the solar neutrino data, including the day- 
night asymmetry, should depend rather weakly on this parameter, this dependence may 
not be negligible. For example, the current 3a limit sin^ ^13 < 0.054 [9] means that the 
day-night difference of the solar neutrino signal can be suppressed by up to 15% compared 
to the 2-flavour (i.e. ^13 = 0) case. If one fixes Amg^ = 2 x 10^^ eV^ (which is the 
best-fit value coming from the analysis of the atmospheric neutrino oscillations performed 
by the Super-K collaboration [32]), then the 3a upper bound is sin^ ^13 < 0.066 [9], and 
the day-night difference of the solar neutrino signal can be suppressed by as much as 20% 
compared to the 2-flavour case. 

While the accuracy with which the solar neutrino day-night effect is measured by the 
current experiments is insufficient for probing the value of ^13, future very large water 
Cherenkov detectors, such as UNO or Hyper-Kamiokande, may be able to provide a sig- 
nificant information on it. Depending on the measured value of A^o, the current upper 
bound on ^13 may be improved, or even a lower bound on this mixing parameter may be 
obtained. An advantage of the day-night asymmetry as opposed to the day-time signal is 
that independent of the overall normalization of the solar neutrino flux, which is 
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currently not known with sufficient accuracy. 

It should be noted that the Earth matter effects on the oscillations of the supernova 
neutrinos inside the Earth are governed by the same quantity, — p!^^ , which determines 
the Earth matter effects on the solar neutrino oscillations (see, e.g., [33]). Therefore our 
results, Eqs. (|3.5D and ( |3.(]| ), can also be used for studying the supernova neutrino oscil- 
lations inside the Earth. Note, however, that the typical energies of supernova neutrinos 
{E ~ 30 MeV) are higher than those of solar neutrinos, and so the expected accuracy of 
our approximation for supernova neutrinos is ~ (10 - 15)%. 
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Note added. As we completed our paper, two new papers appeared [34, 35], in which 
analytic formulas for the Earth matter effect on neutrino oscillations, similar to ours, were 
derived. The first [34] is based on adiabatic perturbation theory, while Ref. [35] employs 
ordinary perturbation theory, similar to that used in the present paper. Note, however, 
that Refs. [34, 35] adopt the two-flavour approximation, while our analysis is performed in 
the three-flavour framework. 



A. Derivation of Eq. (IgTel) 



We shall derive here expression ( |3.£| ) for the Earth regeneration factor — P^^ ■ In 
contrast to perturbation theory in V developed in sec. ^ we employ here perturbation 
theory in the small parameter V^/2(5, so that the obtained results will be valid for arbitrary 
values oiVL. 

Consider the neutrino evolution equation in the mass-eigenstate basis, i{d / dt)i'mass = 



mass'^mass 



explicitly as 



Here Ur, 



Wv, and the Hamiltonian H„ 



WHU can be written 



Hr, 



.2 V 



Sl2 Cl2 Cjg 

V ci2 si3 ci3 e^-^cp V 
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*12 ''13 

si2 si3 ci3 e*^cp y 
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Since V <^ 25 <^ 2A and S13 ^ 1, one can neglect the (1-3) and (2-3) elements of Hmass 
compared to the (3-3) element. This means that the evolution of the third mass eigenstate 
f3 essentially decouples from that of the first two eigenstates and, in addition, that matter 
effects on the third mass eigenstate are negligible. 

Diagonalization of the (z^i, z/2)-subsector of Hmass is performed through the rotation 
by the angle 0'{t) which is determined according to 



cos 26' {t) 



8 -cos 2ei2cl^V{t)/2 



sm2e'{t) 



sm2ei2cl^V{t)/2 



(A2) 



uj{t) ' " uj{t) 

with uj{t) given in Eq. ( 2.17] ). Note that 6' = §12 — Ou, where O12 was defined in Eq. pl^ ). 
The smallness of the potential V compared to 26 means that the neutrino matter eigenstates 
almost coincide with mass eigenstates, i.e. 9'{t) <^ 1. We shall now employ perturbation 
theory in sin 6', which is essentially the same as perturbation theory in V/26. 
The decoupling of 1^3 allows one to write 

I -cos 26' {t) s\n26'{t) \ 



u:{t) 



\ 



sm26'{t) cos26'{t) 

2A/io{t) J 



(A3) 



Here A = A — {V/2 + S)/2, and we have redefined the common phase of the neutrino 
states so as to make the Hamiltonian of the (1^1, 1^2) sector traceless and used Eq. ( [A2D . 
The Hamiltonian Hmass in Eq- (A3) can be decomposed into the term of zeroth order in 



s' = sin 6' and term that contains first and higher orders in s' according to 



f -uj{t) \ 



Hr, 



uj{t) 

2A y 



/ s'(t)2 s'{t)c'{t) 0\ 



+ 2L0{t) 



s'{t)c'{t) -s'{tf 






0/ 



ml 



(A4) 



To leading order in s'(t), only the (1-2) and (2-1) elements of H„ii are different from zero. 
The evolution matrix to first order in s' can be found from a formula which coincides with 
Eq. (3J) with Sq and Hi replaced by [Smass)o and i^mi, respectively. The zeroth-order 
evolution matrix in the mass eigenstates basis is {Smass)o = diag(e*''^, e"**^, /), and to first 



order in s' we find 



-iCe'* \ 



-iC*e-^'t> 




-i<j> 



where 



Loit')dt' 



to 



C = 2 I 'a;(0s(0e"''^*o"^"^''"dt'. 
'to 



(A5) 



(A6) 
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The parameter / is not relevant to our discussion. The amplitude A®(z^2 — ^ i^e) is obtained 
as 



ci3(.-ici2Ce"^ + 3126-''^) . (A7) 



A®{u2^jye) = {l 0)USmass 1 

\v 

This gives 

P2e - P2e = 4 2^12 lm{e^''*'C) , (A8) 

which immediately leads to Eq. ( |3.6| ). 

The two-flavour limit of Eq. (|3.6| ) was derived in Ref. [35] by dividing the interval [0, L\ 
into small intervals of constant matter densities. 
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